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Abstract 



A sticky spheres model to describe slow dynamics of a non-equilibrium sys- 
tem is proposed. The dynamical slowing down is due to the presence of 
entropy barriers. An exact steady state analysis of the representative mean 
field equations, in the case when the clusters are chosen with the same a 
priori probability demonstrates a non-equilibrium phase transition from an 
exponential cluster size distribution to a powerlaw. 



05.40.+j, 64.60.-i, 68.45.Gd, 75.10.Nr 



Typeset using REVTgX 



1 



I. INTRODUCTION 



When a macroscopic system in equilibrium at high temperature is quenched rapidly to 
low temperature, either or both of the following can happen. The system may get thermally 
arrested in a metastable (local energy minimum) state much faster than it could equilibrate 
at that temperature and hence, its subsequent dynamical evolution becomes slow, for details 
see [1]. On the other hand, the system may still have thermal freedom to sample a large 
number of equal or almost equal energy states, upon a temperature quench, so that its 
dynamical evolution again becomes slow, for example see [2]. In other words, the system 
would remain trapped for a long time due to the presence of energy and/or entropy barriers. 
As a result, the relaxation of the system to its equilibrium could become anomalously slow. 
It is often history dependent, usually referred to as 'aging', and could become progressively 
slower with time. Glasses [3], obtained by the rapid quenching of liquids, provide simple 
examples of aging systems which evolve slowly forever towards their putative equilibrium 
states; granular systems whose density compaction is logarithmically slow in response to 
mechanical tapping [4] and react ion- diffusion systems [5], provide other recent examples. 
Quite often one finds that these systems develop a certain degree of spatial disorder as well. 
Experimental evidence for such a scenario has recently been reported in the literature [6]. 
An interesting problem in this context is to see whether simple local dynamical rules could 
be devised so as to capture the essential features of the non-equilibrium slow dynamics. 
In particular, it would be of interest to devise dynamical rules that could lead to slow 
logarithmic growth of length scales often found in several systems. To this end, we propose 
in this paper a sticky sphere model to describe the slow dynamics of a non-equilibrium 
system. 

The model consists of hard spheres placed randomly on a regular lattice. The energy of 
the system is defined in such a way that nearest neighbor contacts between the spheres are 
energetically favored, hence the name 'sticky' spheres. An appropriate length scale for this 
system is the mean cluster size. We present numerical evidence to show that this quantity 
grows logarithmically with time at zero temperature. However, for nonzero temperatures, 
it saturates asymptotically to a stationary value. The model and the simulation details are 
discussed in section II. A general mean field formulation of this model is presented in section 
III. An exact steady state analysis of this mean field model for the special case when the 
clusters are chosen with the same a priori probability, described in section IV, shows a phase 
transition from an exponential to a power law cluster size distribution. A brief summary of 
the results is presented in section V. 

II. MODEL 

Consider a regular one dimensional lattice of size, M+N, consisting of N sites unoccupied 
and M sites occupied by hard spheres of size equal to the lattice spacing. Therefore, the 
spheres on nearest neighbor sites touch each other. We assume periodic boundary conditions. 

Let us define the 'energy' of the system, E(t), at time t as the negative of the total 
number of nearest neighbor contacts: 
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where Ck(t) is the total number of k — mers (i.e., clusters consisting of k spheres touching 
each other at time t). We assume that the number of spheres in the system is conserved: 
Tk = ikc k (t) = M. The lowest energy state of the system corresponds to having a single 
M — mer with energy, Eq = — (M — 1), and may henceforth be called the 'ground state' 
of the system. On the other hand, the highest possible energy realizable for the system 
depends on the values of both M and N. 

For given M and N, we can always have a configuration of spheres with a maximum of 
(M — N) nearest neighbor contacts, with M > N. This implies that the maximum energy 
the system can have is given by, E max (M, N) = — (M — N). However, when the system 
consists of only monomers, which can be realized when M < N, the energy is zero. Thus 
we have, 



From Eq. (1), it follows that the energy per particle, e(t) = — 1 + [C(t)/M], where C(t) = 
Sitli c k(t), is the average number of clusters at time t. Since M/C(t) is just the mean cluster 
size, A(i), we have e(t) = — 1 + (1/A(£)), or equivalently, A(t) = 1/[1 + e(t)]. Thus, we have 
a 'sticky' sphere system in which nearest neighbor contacts are energetically favored. 

We start from an initial (t = 0) configuration of the sticky spheres placed on a one 
dimensional lattice segment in such a way that the system is in the highest possible energy 
state for given M and N. At any instant of time t, we choose a k — mer with a pre assigned 
probability, p k . Usually we take p k as k/M, implying thereby that we choose a sphere at 
random with the a priori probability, 1/M. If we have chosen a monomer (k — 1), then it 
can hop either to its left or to its right with equal probability. On the other hand, if we 
have chosen a k — mer with k > 1, then we choose one of its edge spheres (or equivalently, 
edge particles) at random with equal probability. We call it the 'active' particle. We note 
that there is at least one empty site available for the active particle to hop. Consider the 
situation where we have chosen the leftmost sphere of the k — merik > 1) as the active 
particle. This particle can hop to the left. Let there be an I — mer(l > 1) located to the left 
of the active particle such that there are n empty sites in between them. If n — 1, we simply 
move the particle into the available empty site because it does not cost energy. At the end 
of this move, we have an (/ + 1) — mer and a (k — 1) — mer separated by one empty site. If 
n > 1, then we have two possibilities for the particle to hop, as illustrated in fig. 1(a), and 
described below. 

(i) Hopping to the nearest neighbor empty site: If we move the active particle to the 
nearest empty site, then we would be creating a monomer in the system. This process 
would therefore cost one unit of energy. Hence, in order to take care of this energy cost, 
we move it to the nearest empty site with probability e _/3 , where f3 is the inverse of the 
temperature. At the end of this move, we will have a monomer located in between an 
I — mer and a (k — 1) — mer. 

(ii) Hopping to the farthest empty site: If the above move is not accepted, then we move 
the active particle to the farthest empty site so that it sticks to the right edge of the I — mer. 
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The energy of the system (or equivalently, the number of nearest neighbor contacts in the 
system) does not change. At the end of this move, we will have an (/ + 1) — mer and a 
(k — 1) — mer, with n empty sites in between them. 

We have simulated the above process for the case M = N so that the dynamics will 
cover the full range of energy (e = E/M) from to —1. The a priori probability, pt, for 
choosing a k — mer is taken to be proportional to k in the simulation. We have presented in 
Fig. 2 the mean cluster size, A(t), as a function of ln(t) obtained by averaging the data over 
50 independent runs for a system of size N = 16384 and temperatures j3 = 2,4,5,6,8, 10 
and oo. We observe that A(t) saturates asymptotically for temperatures T > 0, whereas it 
continues to grow logarithmically at T = 0. 

We note that the logarithmically slow dynamics at zero temperature is purely due to 
entropy barriers because monomer creation is not possible at this temperature; the system 
evolves only by the process of hopping to the farthest neighbor empty site, which does not 
cost energy. In this sense, our model belongs to the same class of mean field models as that 
of Ritort [2]. In fact, we could anticipate this on heuristic grounds: 

The system will necessarily have to be in the configurational state consisting of a 
monomer and an (N — 1) — mer before it might be able to reach the ground state by 
choosing the monomer with probability 1/N. This is a rare event because the larger cluster 
will always loose a particle with more probability. Precisely the same situation prevails 
[7] in the Ritort's model as well. Hence, we have also shown in Fig. 2 the growth of A(t) 
obtained from Godreche-Luck (GL) mean field formalism [8] of the Ritort's model as con- 
tinuous lines. We observe that our simulation data agree more or less with those of the GL 
values for asymptotic times t > tql, where we have schematically shown the temperature 
dependence of tql in the inset of Fig. 2. Clearly, tql — > oo as j3 — > oo, and the simulation 
data fall on a line parellel to but below the GL line. The sticky sphere system, therefore, 
admits of a mean field description that incorporates the GL formalism at appropriate limits. 

III. MEAN FIELD FORMULATION OF A ONE DIMENSIONAL STICKY 

SPHERE SYSTEM 

The hopping of a single particle to its nearest /farthest neighbour empty site can be in- 
corporated easily in a mean field description by considering the dual representation obtained 
by replacing particles by holes and holes by particles, k — mersik > 1) of the sticky sphere 
system S (Fig. la) correspond to empty intervals of length k in its dual representation S* 
(Fig. lb), and vice versa. The energy of the system is still given by Eq.(l) except that Ck(t) 
now stands for the number of successive empty sites of length k in S*. 

Consider a k — mer, K^, in S having the empty intervals I m and I n to its left and right 
respectively (Fig. la). This corresponds to the empty interval 1% between an m — mer, K* m , 
and an n — mer, K*, in S*. Let P be the rightmost particle of K^. The hopping of P 
to its right nearest neighbour site, Q, in S corresponds to the dissociation of the leftmost 
particle of K* in S*. On the other hand, hopping of P to the farthest neighbour site, R in S 
corresponds to the cluster K* moving as a whole to the left by one lattice unit in S*. Thus, 
the nearest /farthest neighbour hopping of a particle in S corresponds to (single particle) 
dissociation/movement of a cluster in S*. 
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In general, these processes may occur with probabilities q\ and q 2 respectively. For 
convenience, we may rescale the time so as to have these events (namely, single-particle 
dissociation/movement of a cluster) occur with the rates unity and uj = q 2 /qi respectively. 
Spatial correlation in the system may be ignored by treating the k — mers(k > 1) in S* 
as point masses occupying single lattice sites only. This leads to a simplified mean field 
description of the system in terms of a distribution of masses on N lattice sites. We study 
the stochastic evolution of the system in the thermodynamic limit, M, N — > oo with the 
mass density, p = M/N, remaining finite. 

Let fk(t) be the probability that a site will have mass k at time t. By definition, 
J2t=ofk(t) — 1 an d J2kLo^fk('t) — P- Let p k be the a priori probability for choosing a 
k — cluster and, if chosen, let dk be the a priori probability for moving it by one lattice unit. 
The evolution equation for fk(t) can now be written as 

dfk>2(t) 



dt 



= n(t)f k -i(t) - [tt(0 + Xp(t)p k ]fk(t) + Xp(t)p k+ ifk+i(t) 



- UJ {[ Pk d k + A(t)} f k {t) - E k n= lPnd n f n (t)fk-n(t)} , (3) 

where, 

oo oo 
7T(t) = Y,Pnfn(t); A(t) = J2Pndnfn(t); 
n=l n=l 

oo 

X p (t) = (1 - e-P) 8 (t) + e-' 3 ; s(t) = ^ f n (t) (4) 



n=l 



This equation consists of two parts, one corresponding to the single particle dissociation 
and the other to the cluster moving by a lattice unit as a whole. Each part has both the 
gain and the loss terms. 

In the case of single particle dissociation, there are two gain terms. The first one corre- 
sponds to the event of a dissociated particle sticking to a (k — 1) — cluster. The second one 
corresponds to a particle dissociating from a (k + 1) — cluster, taking care to account for 
the energy cost, e _/3 , involved in the event of its becoming a monomer. Similarly, the first 
of the loss terms corresponds to a dissociated particle sticking to a k — cluster. The second 
one corresponds to a particle dissociating from a k — cluster, taking care to account for the 
energy cost, e _/3 , in the event of its becoming a monomer. The probability of choosing a 
k — cluster, p k {k > 1), has been introduced appropriately. 

In the case of a cluster moving by one lattice unit as a whole, the gain term corresponds 
to an n — cluster(l < n < k) coming to stick to a (k — n) — cluster . The event of a k — cluster 
moving out as well as that of a cluster coming in to stick to a k — cluster constitute the loss 
terms. The probability of moving a cluster, p n d n {n > 1), has been introduced appropriately. 

Similarly, the master equations satisfied by the fractions, f (t) and fi(t), can be written 
as follows: 

dh{t) 



dt 



Hp(t)fo(t) ~ br(t) +pi]/i(t) + Mt)p 2 f 2 (t) 
- u {[pxd! + A(*)]/i(t) - Pidi/i(t)/„(*)} , (5) 

df (t) 



dt 



= -fip(t)f (t) +us(t)A(t), (6) 
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where, Hp(t) = p±(l — e~^)fi(t) + 7r(t)e _/3 . In this model, the parameters 00 and gPs are all 
assumed to be temperature independent. 

The mean field equations, obtained by ignoring the spatial extensions of k — mers(k > 1) 
in S*, provide the simplest representation of the nearest /farthest neighbour single particle 
hopping of a sticky sphere system S. Yet, we can not assume a priori that they describe the 
asymptotic dynamical behaviour of S. Because, the probability, pk, of choosing a k — mer 
actually stands for the probability of choosing the empty interval bounded on one side by 
the k — mer of interes. It is also important to note that the presence or absence of the 
aggregation term, F k = Y,n=iPnd n fnfk-n, in Eq.3 corresponds to the specific monomer 
dynamics implemented in S, viz., whether they jump to their farthest or to their nearest 
neighbour sites respectively However, in the case when the clusters are chosen with equal 
a priori probability, empty intervals are also chosen with the same a priori probability 
(say, pk = 1); hence, the mean field equations could provide an adequate description of S. 
Moreover, it turns out that an exact steady state analysis of these equations can be carried 
out in this case. 



IV. STEADY STATE ANALYSIS 

Here, we consider the case pj = dj = 1, for which an exact steady state analysis can be 
carried out. It is clear from Eq.(3-6) that the generating function, Qp(z,t) = Y^k=i zk fk(t), 
satisfies the following equation: 



dQp(z,t) _ 2 ( , 
dt ~ Qf}{) 



a(z) + 



Qb + c{z) 



(7) 



where, 



A/3 sz 



a(z) = 2s H h 
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b( Z ) : 

c{z) = \p(z - 1) 



(8) 
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In order to study the steady state behavior of the system, we set dQp{z,t)/dt = and 
choose the root of the resulting quadratic equation so that Qp{z = 0) = is ensured: 



2Q p {z) = [a(z) + -}- 

z 



\ 



a(z) + - 

z 



-Ac(z). 



(9) 



Simplifying the algebra, we can show that 



a(z) + 



Ac{z) = 



s(z-l) 



idZ 



(z - zi)(z - z 2 ), 



(10) 



where, the roots, z^ 2 are given by 
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^ + l-Tj[l + 2wT 2y/u 2 + uj] ; r = e _/? . (11) 



Hence, we have the generating function, 

2us + s + \p \p sz s(l-z) r- —— — 

QM = %o 2^z ~ 2^ + -%Jz~^ Z ~ Z1){Z ~ Z2) (12) 

The value of s is fixed by the conservation of particle density, p: 

l = 't[ Xp ~ 8 ( 1 + V / ( 1 -^)( 1 "^))] (13) 

For a given u, it is clear that the value of z±, being always less than z 2 , should not be less 
than unity for p to be real. As p increases, the steady state value for the number of clusters, 
s, increases, thereby reducing the values of z 1}2 - Hence, we have the condition, 




The equality sign defines the critical value s c at which the root z — 1, and hence the critical 
density, 

p c = M / ( P } , ; pI = Ji + ± - 1 (15) 

The number of clusters will not increase beyond s c for p > p c . It is of interest to consider 
the question of how this inequality influences the cluster size distribution. To this end, we 
consider the following contour integral, 

A = (16) 

The contour is chosen suitably so that only the portion of the contour above and below 
the branch cut z — Z\ contributes to the integral. The number of k — mers, f^, has the 
asymptotic exponential form, [l/zi) k for p < p c whereas it has a power law form, k~ 5 ^ 2 
for p = p c ; as the density is increased beyond p c , in addition to the power law decay, the 
distribution develops a delta function peak corresponding to an 'infinite' aggregate. 

However, at zero temperature, z\ < 1 for all nonzero values of u; therefore, the above 
steady state analysis breaks down. In fact, the condition expressed by Eq. (14) can be 
rewritten as 

T>Ta r = ^-f^)i (17) 

For a given a;, the value of r c increases as we increase the particle density, until it becomes 
equal the given temperature; beyond this, the steady state analysis breaks down. In other 
words, the steady state phase transition from the 'exponential' regime to the 'aggregating' 
regime is observable only in a limited range of temperature decided by u> and p. The infinite 
temperature version of a related model has been discussed by Majumdar, Krishnamurthy 
and Barma [9]. 
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V. SUMMARY 



In this paper, we have presented a generic sticky sphere model for describing the non- 
equilibrium behavior of a system fast quenched to a low temperature. The evolution of 
the system is based on a local dynamical rule - the nearest/farthest neighbor hopping of 
a randomly chosen particle. The mean cluster size, defining a length scale for the system, 
asymptotically saturates to a stationary value at nonzero temperatures, whereas it grows 
logarithmically with time at zero temperature. We have presented a general mean field 
formulation of this model and solved it exactly for the case when the clusters are chosen 
and moved with the same a priori probability. We have shown that the steady state cluster 
size distribution undergoes a phase transition (in appropriate temperature range) from an 
exponential form to a power law with an additional delta function peak corresponding to an 
'infinite' cluster. 

One of the authors (S.L.N) thanks Y. S. Mayya and Amitabh Joshi for fruitful discussions. 
*glass@apsara.barc.ernet.in; 
+ kpn@igcar.ernet.in 
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FIGURES 




FIG. 1. (a)Sticky sphere model: the dynamical moves available for a chosen particle; (b) its 
dual representation: the corresponding dynamical moves, r = e _/3 . 
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FIG. 2. Mean Cluster size for N = 16384. The time t is measured in units of 1/N. Inverse 
of temperature j3 = oo, 10,8,6,5,4 and 2, from top to bottom. Open circles represent simulation 
data obtained as 50 runs averages; Continuous lines have been obtained from Godreche and Luck 
mean field formalism[4] of the Ritort's model. Inset: Temperature dependence of tql, the times 
beyond which simulation more or less agrees with GL. 
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